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We develop the theory of AbeHan functions defined using a tetragonal curve of 
genus six, discussing in detail the cyclic curve = + X4x'^ + X^x^ + X2X^ + 
Xix + Aq. We construct Abelian functions using the multivariate cr-function 
associated to the curve, generalising the theory of the Weierstrass p-function. 
We demonstrate that such functions can give a solution to the KP-equation, 
outlining how a general class of solutions could be generated using a wider class 
of curves. We also present the associated partial differential equations satisfied 
by the functions, the solution of the Jacobi Inversion Problem, a power series 
expansion for a{u) and a new addition formula. 

1 Introduction 

Recent times have seen a revival of interest in the theory of Abelian (multiply 
periodic) functions associated with algebraic curves. This topic can be dated 
back to the Weierstrass theory of elliptic functions, which we use as a model. 
Let a{u) and p{u) be the standard Weierstrass functions (see for example [21]). 
The p-function can be used to parametrise an elliptic curve, = Ax^ — g2X — g^, 
and satisfies the following well-known formulae. 



(1) 




(2) 
(3) 



The (T-function satisfied a power series expansion 
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5253"" + ■ • ■ (4) 



and a two term addition formula 

a{u + v)a{u — v) 



p{u) - p{v). 



(5) 
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Taking logarithmic derivatives of this will give the standard addition formula for 
p(u). This paper will generalise equations (l)-(5) for a previously unconsidered 
class of functions.. 

The study of Abclian functions associated with the simplest hyperelliptic 
curves (those of genus 2) goes back to the start of the 20th century. Klein's 
generalisation of the Weierstrass theory is described in Baker's classic texts, [1] 
and [2], while [5] gives a more recent study of the general hyperelliptic case. 
Further generalisation has been structured by considering, (with the notation 
of [8]), classes of (n, s)-curves. These are curves with equation 

- - ^ fJ.[ns-an-f3s]X"y'^ fJ^j COUStautS, (6) 

a, (3 

where a, /3 G Z with a e (0, s — 1), /? G (0, n — 1) and an + f3s < ns. The cyclic 
subset of such a class of curves is generated by setting [3 — 0. We suppose that 
(n, s) are coprime, in which case the curves have genus g = ^(n — l)(s — 1), and 
a unique branch point oo at infinity. 

In the last few years a good deal of progress has been made on the theory 
of Abclian functions associated with trigonal curves, (those with n = 3). The 
cr- function realisation of these functions was developed first in [7] and [12], with 
the two canonical cases studied in detail in [10] and [3]. 

In this paper we consider the next logical class, and work with a tetrag- 
onal curve, (one with n = 4). We have started by looking at the curves of 
lowest genus, and simplified by considering the cyclic subclass. We construct 
the multivariate cr-function associated with this curve, and use it to define and 
analyse classes of Abelian functions, generalising the theory of the Weierstrass 
p-function. A key component of our work was the construction of a scries ex- 
pansion for the (T-function. This technique was first developed for the trigonal 
case in [4], however the computation involved for the present expansion is sig- 
nificantly greater. The latter computations were performed in parallel with the 
use of the Distributed Maple software, (sec [19] and [20]). 

The applications of Abelian functions to integrable systems and soliton the- 
ory have been the topic of research for some time, (see for example [16] and 
[13]). It is well known that the elliptic p-function could be used to construct 
a solution to the KdV-equation. Similar solutions to non-linear equations have 
been derived from higher genus curves, for example in [12], where the function 
p33 associated with the (3,4)-curve was shown to be a solution of the Boussinesq 
equation. This has suggested a more general link between such functions and 
the integrable KP hierarchy. We have demonstrated how the Abelian functions 
we define can give a solution to the KP-equation, outlining how similar solutions 
will also be found from any (4, s)-curvc. 

This paper is organised as follows. We give the basic properties of the curve 
we consider in Section 2, including explicit constructions of the differentials on 
the curve, and a set of weights that render the key equations homogeneous. Then 
in Sections 3 and 4 we define the cr-function and Abelian functions associated 
with this curve. Section 5 discusses a key theorem satisfied by the p-functions, 
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which we use to give a solution to the Jacobi Inversion Problem. In Section 6 we 
derive some properties of the cr-function, including the series expansion, while 
in Section 7 we use this to generate relations between the Abelian functions. 
Section 8 demonstrates how solutions to the KP-equation can be constructed 
from Abelian functions. Finally, in Section 9 we give the derivation of a two- 
term addition formula. 

2 The purely tetragonal curves 

We will investigate Abelian functions associated with a tetragonal curve. The 
simplest general tetragonal curve is, in the notation of the (n, s)-curves, a (4,5)- 
curve. This would be given by g(x,y) = where 

9{x, y) = + {^^lx + fj.5)y^ + {fj-2x^ + fj-ex + mio)?/^ 

4- {^i3X^ + + mix + i-H5)y 

— [x^ + fj./ix'^ + i-igx^ + fii2x'^ + Hiqx + fj.2o) il^j constants). 

In this paper we further simplify by considering the cyclic subclass of this family. 
These are the curves C, given by 

C : f{x,y) = {Xj constants) 

where f{x, y) = / - {x^ + X^x'^ + X^^x^ + Azx^ + Xix + Aq) . (7) 

The curve C has genus g = 6, the unique branch point oo at infinity and is 
referred to as the purely tetragonal^ or strictly tetragonal curve. It contains an 
extra level of symmetry, demonstrated by the fact that it is invariant under 

[C]: {x,y)^{xXy). (8) 

where C is a 4th root of unity. 

For any (n, s)-curve we can define a set of weights for the variables of the 
theory, including the curve constants, which render equation homogeneous with 
respect to the weights. To find these weights consider the mapping x ^ i"^X 
acting on all elements in the curve equation. Then define the weights as the 
constants ai that render the new equation homogeneous with respect to t. The 
weights of x,y will then be determined up to a constant by nay = sa^- To 
keep with convention, we let ax = —n and ay = — s so that they are the largest 
negative integers satisfying this condition. The weights of the curve constants 
can then be chosen to make the remainder of the equation homogeneous. 

Definition 2.1. For the cyclic (4,5)-case we have 





X 


y 


A4 


A3 


A2 


Ai 


Ao 


Weight 


-4 


-5 


-4 


-8 


-12 


-16 


-20 



while in the general (4-,5)-case the weights of the curve constants are given by 
their subscripts. We refer to these as the Sato Weights 
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As we precede through the paper we can use the approach of this mapping to 
conclude that other elements in our theory must have definite weight, and all 
the equations presented here will be homogeneous with respect to these weights. 
Next we construct the standard basis of holomorphic differentials upon C . 



du = (dwi, . . . , rfug), duAx, y) = dx, 

4yJ 



Where 9i{x,y) = l 52(x,y)=x, g,{x,y)^ y 

g4{x,y)=x^, g5{x,y}^xy, gG(x,y) = y^. 

Denote points in by u for example, and their coordinates by (ui, M2, . . . , ug). 
We know from the general theory, that any point u E can be expressed as 



6 /.P, 

(ui,M2,U3,U4,ti5,M6) = ^ / du, 

1=1 



where the Pi are six variable points upon C . Let A denote the lattice generated 
by the integrals of the basis of holomorphic differentials along any closed paths 
in C . Then the manifold C^/A is the Jacobian variety of C, denoted by J. Let 
K be the map of modulo A over C: 

K : ^ CVA = J. 

Therefore A = k"^((0, . . . , 0)). Next, for k = 1,2,... define 21, the Abel map 
from the /cth symmetric product Sym'^(C) to J. 

21 : Sym'=(C) ^ J 

(Pi,...,Pfc) ( / du+-+ I du] (mod A). (10) 




where the Pi are points upon C. Denote the image of the fcth Abel map by 
t^I*^], and let 

[-1](mi, ...,uq) = (-ui, . . . , -lie). 
Define the kth standard theta subset (often referred to as the fcth strata) by 

When fc = 1 the Abel map gives an embedding of the curve C, upon which we 
define ^ as the local parameter at the origin, 2li(oo). 

£, = .T^" = (11) 

We can then express the basis (9) with ^ and integrate to give 



(12) 
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The higher order terms will contain the curve parameters A = {Aq, . . . , A4}. 
Similarly to [4] and [12], we could rewrite this using as the local paremeter. 

Note that such calculations can be performed similarly for any (n, s)-curve, 
and that since each element of du is homogeneous in Sato weight we can con- 
clude that the m have definite Sato weight. Since the weight of ^ must be +1 
from equation (11) we can define the weights of u uniquely as below. 

Definition 2.2. In the (4-,5)-case we assign the following weights to u. 





Ui 


U2 


"3 


U4 


U5 


Ui 


Weight 


+11 


+7 


+6 


+3 


+2 


+1 



Remark 2.3. The weights of the variables coincide with the order of their zero 
at 00. They can also be calculated using the Weierstrass gap sequence, where 
the weights of ui, . . . ,uq are the gap numbers, and the weights of x and y are 
the negative of the first two non-gap numbers. 

Definition 2.4. Let {x,y) and {z,w) be two variable points upon C. Then the 
2-form Q(^{x,y), {z,w)) on C x C is a fundamental differential of the second 
kind if 

1. It is symmetric: n[{x,y), {z,w)) = n[{z,w), (x,y)) . 

2. The only pole of second order is along the diagonal of C x C (where x = z). 

3. It can be expanded in a power series as 

n{{x,y),{z,w))^(^j^^ + 0{l)^d^d£,' ( as {x,y) iz,w)), 

where ^ and ^' are the local coordinates of {x,y) and {z,w). 

We will construct Klein's explicit realisation of this in Proposition 2.6 below. 

First introduce dr, the basis of mcromorphic differentials which have their 
only pole at 00. These are determined modulo the space spanned by the du 
and can be expressed as 

h ' (x Ui 

dr = {dri, . . . ,drQ), where drj{x,y) = ^ ^ — dx. (13) 

An explicit basis is constructed later, in order to satisfy Proposition 2.6. 
Definition 2.5. Define the following meromorphic function on C x C as 

where [ ]w means that we remove any terms which have negative powers with 
respect to w. 



yji-k+l 
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Proposition 2.6. The fundamental differential of the second kind can he ex- 
pressed as 

^{{.x, y), {z, w)) = R{{x, y), {z, w))dxdz, 

where 

The polynomials hj{x,y) need to be chosen so that D, is symmetric. This will 
lead to a realisation of fl in the form 

nil ^ / ^^ F {{x , y) , (z , w)) dxdz 

Proof: The essential part of the proof is the same as in the lower genus 
cases (see [10] for example). In this case we explicitly determine the basis of 
meromorphic differentials, (13), to be given with 

hi ^ -y'^{%x^\i + llx^ + ^x\z + 2A2), /i2 = -y'^{h + 4xA4 + 7x'^), 
h^ = — 2a;j/(A3 + 3x^ + 2xX4^, /i4 = —3xy^, /15 = ~2x^y, /ig — —x^. 

The polynomial F in the realisation (14) is found to be, 

F{{x, y), {z, w)) = Ay^w^ + (3a;z* + z^Ag + z^x^ + 2X2Z^ + SxAgz^ 
+ Az^xXi + 4Ao + Xix + 2\2xz + 3Aiz)y^ + (2Aiz + A\2xz + 4Ao 
+ 2Aia; + 4x^A4Z^ + 2A3X^z + 2x^z^ + 2z^x^ + 2xX3z'^)wy 
+ {X^x^ + 4Ao + 3Aia: + 2A2a:^ + Aiz + x^z'^ + 3x^z + 2A2XZ 
+ 3X3X^z + AXiX^z)w'^. (15) 



3 Defining the <T-function 

In this section we describe the multivariate cr-function associated with C, from 
which all Abelian functions associated with C can be defined. This can be re- 
garded as a generalisation of the Weierstrass cr-function, with the main difference 
that there are now g = 6 variables. 

(T = a{u) = a{ui,U2,U3,U4,U5,ue). 

First we choose a basis of cycles (closed paths) upon the surface defined by C. 
We denote them 

ai,f3j, l<i,j<6, 
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and ensure they have intersection numbers 



f3^ ■ Pj = 0, 



ai ■ /3. 



1 if i = j 
if i^j 



This allows us to define the following period matrices: 



du. 



dri 



We combine these into 



kj=l,...,6 
k,i=l,...,6 

M = 



I'll. 



drt 



kj=l,...,6 
k.l=l,...,6 



v' v" 



which we know from classical results to satisfy 
M 





T 


'J UJ"' 


M = 2m 




.^6 





(16) 



This is the generalised Legendre equation, (see [5] pll). We also have that 
{uj')~^ijj" is symmetric with 



Im(((jj') ^^^"^ positive definite. 



(17) 



We now define the multivariate cr-function associated with C. This can be 
constructed using the multivariate 0-fimction, (see for example, [17]). 

Definition 3.1. The Kleinian cr-function associated with C is 

a{u)^a{u-M)^ccyiY,{-\ur]'{uj')-^u^) x e[5]{{uj')-^u'^ \ {uo'y^uj") 

= cexp ( — ^urj' {uo')^^u^^ 



exp 2ml i( 



(m + 5'Y {u')-^u"{m + 5') + [m + 5'Y ((a.')" + 5") 



where c is a constant dependent upon the curve parameters, {Aq, Ai, A2, A3, A4} 
and fixed later (see Remark 6.4)- The matrix S = 



6" 



is the theta function 



characteristic which gives the Riemann constant for C with respect to the base 
point 00 and the period matrix [uj' ,uj"], (see [5] p23-24)- 

In this paper we give some of the most important properties of a{u). How- 
ever, for a more detailed study of the construction and properties of the multi- 
variate cr-function, we refer the reader to [5]. 

Lemma 3.2. Given u G C^, denote by u' and u" the unique elements in MP 
such that 
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Let I represent a point on the period lattice 

e = e'uj' + tu" e A. 

For u, D e and ^ G A, define L[u, v) and x(^) as follows: 

Xi£) - cxp [Tn{2{e'^6" - e"^5') + £'^£")] . 
Then, for alZ 16 G C^,^ G A the function a{u) has the quasi-periodicity property 

cr(tt + £)-xWcxp[L(u+^,£)] -criu). (18) 

Also, for 7 G Sp{12,'Z,) we have 

a{u;-/M) ^ a{u;]\I). (19) 

Proof: The quasi-periodicity property given in equation (18) is a classical 
result, first discussed in [1], that was fundamental to the original definition of 
the multivariate cr-function. Equation (19) is easily seen from the definition of 
(t{u), since 7 corresponds to the choice of basis cycles {ctj, l3j}^^i which were 
used to define M. 

m 

4 Classes of Abelian functions 

Definition 4.1. Let 9JI(m) be a meromorphic function of u E C^. Then 9JI is 
an Abelian function associated with C if 

m{u + u'n^ + Lj"m^) = m{u), 

for all integer vectors n,m €z Z, wherever 9Jl is defined. 

We now define a set of fundamental Abelian functions on J. 

Definition 4.2. Define the 2-index Kleinian p-functions as 

Py(tt) = -— — — logcr(M), i<je{l,...,6}. 

ouiOUj 

A short calculation shows these functions to have have poles of order 2 when 
ct{u) = 0, and no other singularities. We can check (using Lemma 3.2), that 

Pjj(m + £) = P,j(m), V £gA. 

Hence we can conclude these functions to be Abelian. Similar analysis will show 
their derivatives to be Abelian also. 
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Definition 4.3. For n> 2, define n-index Kleinian p-fimctions as 

P,.......„ in) = • . . ^ logaiu), < . . • < z„ G {1, . . . , 6}. 

Remark 4.4. (i) Compare with equation (1) and the elliptic case to see that we 
are defining a generalisation of the Weierstrass p-function, and its derivatives, 
(ii) This notation is compatible with the elliptic case, where we would now de- 
note the Weierstrass p-function as pii(u) and its first derivative p' by pmiu). 
(Hi) The order of the indices is irrelevant. For simplicity we always use ascend- 
ing numerical order. 

(iv) We are usually only referring to one vector of variables u. In these cases, 
for simplicity, we write pij instead of pij{u). 

We find in Section 7 tliat tlie p-funetions are not sufficient to construct a basis of 
the simplest Abelian functions. Hence we also define a generalisation of Baker's 
Q-functions, which we need to extend further than in the lower genus cases. 

Definition 4.5. Define the operator as below. This is now known as Hirota's 
bilinear operator, although it was used much earlier by Baker in [2] for example. 

^ _ _d d_ 

* dui dvi 

It is then simple to check that an alternative, equivalent definition of the 2-index 
Kleinian p-functions is given by 



p,j{u) = --^^--^—:^A,Aja{u)a{v) _ « < i S {1, • • • , 6}. 



1 

2aju) 

We extend this to define n-index Q-functions, for n even. 



= } ^\o ^^Ai.■■■^^^cr{u)c^{v) ii < ... <ine {!,..., 6}. 

2(t(m)^ ■ v=u 

We can show, as above, that these functions are also Abelian. 

Remark 4.6. (i) The subscripts of the p-functions denote differentiation 

d 



du 



"Pil, 42, ■■■,!„ — Pii,i2, 



but this is not the case for the Q -functions. Here the indices refer to which 
Hirota operators were used. 

(ii) If we had applied the definition for n odd, then it would have returned zero. 
(Hi) Note that both the p-functions and the Q-function have poles when a{u) = 
and no other singularities. However, the n-index p-functions had poles of order 
n, while the n-index Q-function all have poles of order 2. 

The 4-index Q-functions were first used by Baker, and in [10] it was shown 
that they could be expressed using the Kleinian p-functions as 

Qi]M = Pt]ki - 2ptjpki - 2pikpje - 2papjk- (20) 



9 



Proposition 4.7. The 6-index Q-functions can be written as 

Qijklmn — 'Pijklmn 2 ^ |^p2j p/^^mn H" 'Pik^jlmn H" 'Pil^jkmn H" ^im^jkln 

+ PinPjkhn] + [pjkpilmn + pjlpikrnji + Pjmpikln + PjnP'ifc/m] 
~l~ [Pklpijmn H~ Pkmpijln H" Pknpijlm m'Pijkn ~t~ ^In^ijkm] 

+ PmnPijkl^ + ^{[pijPklPmn + PiiPfemPln + Pij Pfcrip/m] 

+ [pikPjlPmn + PifcPjmP/n + PifcPjnPlm] + [pilPjkPmn 
+ PiiPjmPfcn + PilPjnPkni] + [pimpjfep/ri + PimPj/Pferi 

+ PimPjnPkl] + [pinPjfcPlm + PinPjlPkm + PiriPj-mPfcij^ (21) 

Proof: Apply Definitions 4.2 and Qdef to reduce the equation to a sum of 
cr-dcrivatives. We find that they aU cancel (Maple is useful here). The structure 
of the sum was prompted by considering the result for the 4-index Q-functions. 

■ 

Clearly, equation (21) will specialise to give a set of simpler formulae, such as 



innnnnn — Pnn7innn '^^Pn7iPnnn7i ~^ ^^Pnn 



5 Expanding the Kleinian formula 

This section is based upon the following Theorem (originally by Klein). It is 
given for a general curve as Theorem 3.4 in [12]. From this theorem we are able 
to solve the Jacobi Inversion Problem, as well as generate relations between the 
p-functions 

Theorem 5.1. Let {Pi, Pq\ £ be an arbitrary set of distinct points on 
C , and (2, w) any point of this set. Then for an arbitrary point {x, y) and base 
point 00 on C we have 

du\ g,[x,y)g,[z,w) ^ (^^3^)2 ' ^ > 

Here gi is the numerator of dui, as given in equation (9), and F is the sym- 
metric function appearing in equation (15) as the numerator of the fundamental 
differential of the second kind. 

We use our explicit calculation of the differentials to construct (22). We 
expand this as one of the Pj. tends to infinity, to obtain a series expansion in 
terms of the local parameter ^, given earlier in equation (11). It follows that 
each coefficient with respect to ^ must be zero for any u £ J and some (z, w) on 
C. This gives us a potentially infinite sequence of equations, starting with the 
five given in Appendix A. The first 14 have been calculated explicitly (using 
Maple), and can be foimd online at [15]. 
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Manipulating these equations 

We follow the approach of the trigonal papers, and manipulate these equations 
using Maple. We first take the resultant of pairs of these equations, (eliminating 
the variable w by choice), to give a new set of equations dependent on z and 
the p-functions. We let Res(a, b) represent the resultant of equations (a) and 
(b). 

These new equations arc considerably longer than those obtained in the 
lower genus cases. We need to combine them to give a polynomial of degree 
(? — 1 = 5 in z. Such a polynomial would have only 5 solutions, but must be 
satisfied for all u (which has 6 variables). Hence all the coefficients must be 
zero, giving us a set of relations between the p-functions. However, we have 
the extra complication, (compared to the trigonal cases), that none of the new 
equations has degree in z equal to g. Therefore, in this case, at least two rounds 
of elimination between the equations will be required. 

We find that Res(34,35) has degree 7 in z, so we rearrange it to give an 
equation for z''. Then since Res(34, 36) and Res(34, 38) have degree 8, we can 
repeatedly substitute for z'^ in both until we are left with two equations of degree 
6 in z. Since these are very long we do not print them here, however they can 
be found online at [15] where we have labelled them (Tl) and (T2). 

We next rearrange (Tl) to give an equation for z^ and repeatedly substitute 
for 2*' and z^ in the remaining equations until they are of degree 5 in z. The 
coefficients of such equations must be zero, giving us relations between the p- 
functions. The smallest such relation has 3695 terms, with the others rising in 
size considerably. Unlike the trigonal cases, these can not be easily separated 
to give expressions for individual p-functions. However, these are implemented 
in the construction of the cr-fimction expansion, (see Section 6). 

Jacob! Inversion Problem 

Recall that the Jacobi Inversion Problem is, given a point m e J, to find the 
preimage of this point under the Abel map (10). 

Theorem 5.2. Suppose we are given {ui,...,Ug} = it £ J. Then we could 
solve the Jacobi Inversion Problem explicitly using the equations derived from 

(34) -(38). 

Proof: Consider either equation (Tl) or (T2) defined in the discussion above. 
This was a polynomial constructed from p-functions and the variable z. This 
equation had degree 6 in z so denote by (zi, . . . , zg) the 6 zeros of the polynomial. 
Next, rearrange (34) to give an equation for w'^. Substitute this into equation 

(35) and multiply all terms by pgg to give the following equation of degree 1 
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with respect to w: 

= w(2:p66p55 - 22^p66 " Zpeep566 + P666^P56 + P36p666 + P66p35 

- zpIq - p36p56 - Peepaee) + z^pmp^h - z^pmp^m + p56^^ 

- 2:p66p266 " P562:p26 + ^:p66p25 + Pl5p66 " P666^^ + P666p462^^ 

+ P6662:p26 - Pl66p66 + P666pl6 " P56p462:^ " P56pl6- (23) 

Wc could substitute each Zi into equation (23) in turn, and solve to find the cor- 
responding Wi- We can therefore identify the set of points {(21, wi), . . . , (zg, Wg)} 
on the curve C which are the Abel preimage of u. 



6 Deriving the properties of cr{u) 

In this section we derive some properties for a(u) and use them to construct 
the series expansion. 

Lemma 6.1. The function a(u) has zeroes of order 1 when u <E G^'^l. Further, 
(t(m) 7^ for all other u. 

Proof: This is a classical result, first discussed in [1], which always holds 
on els-^l. The first part can also be concluded explicitly from the results of 
the previous section. In Theorem 5.2 we discussed how the six roots of the 
polynomial {T2), gave us the Abel preimage of it e J. Now, suppose that u 
is approaching Q^^\ implying one of these roots is approaching infinity. We 
explicitly calculate the denominator of (T2) to be o'(m)^^, using Definition 4.3. 
Therefore, we can conclude that when u descends to we must have cr{u) = 0. 

■ 

Consider u <E 6'^' which, by definition, we can express using points Pk on C as 




Use equations (12) to express u with five local parameters. 

= "1 + ^^i" + • • • + ^^5" + oiei') +■■■+ o{e,'), 

= + 6 + • • • + ^5 + o(ef ) + • • • + o(ef ). (24) 

Suppose we were to take the multivariate resultant of these six equation, elimi- 
nating the parameters f 1 , . . . , ^5 . From the theory of resultants we would be left 
with the unique, (up to a constant), polynomial that must be zero for m G G^^l. 
By Lemma 6.1 this polynomial would be equal to a{u). Further, since it was 
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generated using polynomials of homogeneous weight, we know that (7{u) must 
also have definite weight. 

We can perform this calculation explicitly in the case when A = 0. In this 
case equations (12) simplify to 

and hence equations (24) become finite polynomials. 

We use a multivariate resultant calculation to eliminate ^i, . . . , ^5 and leave 
the polynomial below equal to zero. 

'5'W^4,5 = 8385528^6^ + 3i6"i"5"4 - j^Wg"! - ^ulusU^ - ^U^UaU^ui 

- hl^^ul - Uiul + ^ujul - ^ulul - ^ului + j^ulu2 

+ 5m|u2 + ^U^uj - + 399168 "6^^4 + UAUeulu2 + 

+ 2 U5U3U2 + ^ U5^'"6^"2 + j2 Uq^U2U4 — i W4^U6^W2 + | U4^U6^U5^ 

— i Ui^UfjUc,'^ — ^ U5'*U4tlg"' + U^UqUi — U5^Ul. (25) 

In fact this is just a specific case of the following result for the cr-function. 

Lemma 6.2. Define the canonical limit of the sigma function as the value of 
cr{u) in the case when all the curve constants are zero. In this case the series 
expansion of a{u) about u = (0,0,0,0,0,0) is given by a constant K multiplied 
by the Schur-Weierstrass polynomial generated by (n, s). 

Proof: The result was first stated in [6] , with an alternative proof now avail- 
able in [18]. 

■ 

Note that calculating 51^4^5 as the Schur-Weierstrass polynomial is. computa- 
tionally, far easier than using a multivariate resultant method. 

Corollary 6.3. The function a{u) associated with the (4,5)-curve, is odd with 
respect to u t-^ [—l]u. 

Proof: Fix the matrix M that satisfies (16) and (17). Then the solutions to 
(18) form a one dimensional space over C, (see [3] p456). 

Since we can express (t([— l]u) as a{u + £) for some £ € A, it follows that 
both a{u) and f7([— l]n) satisfy (18). Therefore we have (t([— l]ii) = ka{u), for 
some fc £ C. If we let u = [—i]u then we see have fc^ = ±1, and so a{u) is 
either odd or even with respect to it i-^ [^l]''^- 

We can easily check that 5^4,5 is an odd polynomial from equation (25), 
and therefore by Lemma 6.2 we conclude that o'{u) is an odd function. 

■ 

We now aim to derive a Taylor series expansion for a{u), similar to that of the 
elliptic case in equation (4). This will depend on the variables u — (ui, . . . , uq) 
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and the curve constants A = (A4, . . . , Ao). From Lemma 6.2 we already have, 

(t(m) = K ■ SWa,c, + terms with degree in {A4, . . . Aq} greater than zero, 

for some constant K. Further, we know by Corollary 6.3 that the expansion will 
be odd, and also, since we know that ct{u) has definite weight, we have that the 
expansion is homogeneous in the Sato weights. From equation (6.2) we can see 
that weight is +15. 

Remark 6.4. For simplicity, we fix the constant c in Definition (3.1) to he 
the constant that makes K = 1 in Lemma (6.2). Note that some other authors 
working in this area would define c to he 



where D is the discriminant of the curve C . Note that this constant cancels 
in the definitions of all the Ahelian functions defined from a{u). Hence, any 
relations between such functions are independent of the choice of c. 

We now have enough information to define the following expansion for cr(tt). 

Theorem 6.5. The function cr(u) associated with (7) has an expansion of the 
following form. 

a{u) = Cr{ui,U2,U3, U4, U5, Uq) = Ci5(m) + Cl9(tt) H h Cl5+4n(M) + . . . 

where each Ck is a finite, odd polynomial composed of products of monomials in 
Ui of total weight multiplied by monomials in Xj of total weight 15 — k. 

Proof: The theoretical part of the proof follows [10] and [3]. The rational is 
that although the expansion is homogeneous of weight +15, it will contain both 
Ui (with +ve weight) and Xj (with -ve weight). We hence split up the infinite 
expansion, into finite polynomials whose terms share common weight ratios. 

The first polynomial will be the terms with the lowest weight in Ui. These 
must be the terms that do not vary with A. The indices then increase by four 
since the weights of A decrease by four (see Definition 2.1). 

■ 

By Lemma 6.2 we have C15 = 6*^4.5 as given by equation (25). Using the 
computer algebra package Maple, we calculate the other polynomials succes- 
sively as follows: 

1. Select the terms that could appear in Ck. These arc a finite number of 
monomials formed by entries of u and A with the appropriate weight ratio. 

2. Construct (j{u) as the sum of Ck derived thus far. Then add to this each 
of the possible terms, multiplied by an independent, unidentified constant. 

3. Determine the constants by ensuring (7(m) satisfies known properties of 
the (T-function. 
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• For the first few Ck this was mainly ensuring Lemma 6.1 is satisfied 
(as in the trigonal calculations). 

• For the latter Ck, the coefficients were instead determined by ensur- 
ing a variety of the equations from Lemma 7.3 were satisfied. 

• In addition, those polynomials up to C39 required we ensure <7{u) 
satisfied some of the relations between p-functions obtained from the 
expansion of the Kleinian formula in Section 5. 

The second method is the most computationally efficient (due to the pole 
cancellations), while the third method is extremely difficult. The equa- 
tions in Lemma 7.3 are derived in tandem with the cr-function expansion, 
and so cannot be used for the first few Ck- 

The expansion has been calculated up to and including C59. Appendix B con- 
tains Cig and C23 with the rest of the expansion online at [15]. These latter 
polynomials are extremely large, and represent a significant amount of computa- 
tion. Many of the calculations were run in parallel on a cluster of machines using 
the Distributed Maple package (sec [19]). This expansion is sufficient for any 
explicit calculations. However, it would be ideal to find a recursive construction 
of the expansion generalising the elliptic case, (see for example [11]). 

7 Relations between the Abelian functions 

In the previous section we showed that (j{u) has definite Sato weight, and hence 
so does the Abelian functions defined from it. We can conclude from Definition 
4.3 that 

wt(pii,i2,...,i„) = - [wi(uij -t- wt(uij -t- . . . wi(ui„)] . (26) 
Then use equations (20) and (21) respectively to conclude 

^i{Qi]kl) ^ Vft{pijkl) and Wt((3yfc;™„) = Wt(py7cimn)- (27) 

We now introduce the following definition to classify the Abelian functions 
associated with C by their pole structure. 

Definition 7.1. Define 

r(,/,o(meW)) 

as the vector space of Abelian functions defined upon J, which have poles of 
order at most m, occurring only on the kth standard theta subset, GI*^'. 

Recall that the Abelian functions we define all had poles occurring only when 
(7{u) = 0; which by Lemma 6.1, is when u G G^'^l. Therefore, using Remark 
4.6(iii), we conclude that the n-index p-functions belong to r(j, ©(nGl^l)), 
while the n-index Q-functions all belong to ( J, 0(29^^])) . 
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Theorem 7.2. A basis for T{J, 0(20^^^)) is given by 



CI 

© Cpi5 
® Cp25 
® Cp36 
® Cp56 

© CQ4455 



© Cpii 

© Cpi6 

© Cp26 

© Cp44 

© Cp66 

© CQ3566 



© Cpi2 

© Cp22 

© Cp33 

© Cp45 

© CQssee 

© CQ3556 



© Cpi3 

© Cp23 

© Cp34 

© Cp46 

© 

© CQ2566 



© Cpi4 
© Cp24 
© Cp35 
© Cp55 

© CQ4555 

© CQ2556 



© CQ3456 © CQ2456 © Cg3366 © Cg3445 © CQ2366 

© CQ2445 © CQi466 © CQi556 © CQ2266 © CQ2356 

© CQ2256 © Cg2346 © Cgi455 © Cg2345 © CQ3344 

© CQ2245 © CQ2344 © CQ1266 © CQi356 © CQ1444 

© CQi346 © CQ2236 © CQ2335 © CQi246 © CQ1255 

© CQ1245 © CQ1166 © CQ1244 © CQii56 © CQii46 

© CQ1155 © CQ1145 © CQ1144 © CQii4466- 

Proof: The dimension of the space is 2^ = 2^ = 64 by the Riemann-Roch 
theorem for Abehan varieties. It was shown above that all the selected elements 
do in fact belong to the space. All that remains is to prove their linear inde- 
pendence, which can be done explicitly using Maple. 

■ 

The actual construction of the basis was as follows. We started by including 
all 21 of the pij in the basis, since they were all linearly independent. Then, to 
decide which Qijki to include, wc systematically considered decreasing weights 
in turn, starting at —4 since this is the highest weight of any Q-function. At 
each stage we derived equations to express the Q-functions at that weight using 
the following method (implemented with Maple): 

1. Choose the possible terms at this weight. These are the elements cur- 
rently in the basis with this weight, along with elements in the basis 
of a higher weight (already determined) combined with appropriate A- 
monomials, that balance the weight. 

2. We form a sum of these terms, each multiplied by an undetermined coef- 
ficient. We also include in this sum, the Qijki which are at this weight. 

3. Substitute the Abelian functions for their definitions as cr-derivatives. 

4. Substitute (T(it) for the expansion, truncated at the appropriate point. 

5. Take the numerator of the resulting expression and separate into mono- 
mials in u and A, with coefficients in the unidentified coefficients. 

6. Set all the coefficients to zero, and solve the resulting system of equations. 
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At weights which have more than one Q-function, we often find that one or 
more must be added to the basis, so that the others can be expressed. 

We form these equations at successively lower weights, constructing the basis 
as we proceed. As the weight decreases, we require more of the expansion, 
which is why these were calculated in tandem. Also, as the weight decreases the 
possible number of terms increase, and the computations take more time and 
memory. Upon completing this process we have 63 basis elements. 

We find the final element by considering the 6-index Q-functions. Repeating 
the process, we found that one of the functions at weight —30 is required to 
express the others. 

Sets of differential equations satisfied by the Abelian functions 

We now present a number of differential equations between the Abelian func- 
tions. The number in brackets on the left indicates the weight of the equation. 

Lemma 7.3. Those ^-index Q-functions not in the basis, can be expressed as 
a linear combination of the basis elements. 

(-4) 06666 = -3p55 + 4p46, (-7) 04566 = 2A4P56 + 2p36, 

(-5) O5666 = -2p45, (-7) O5556 = 4A4p56 + 4p36, 

(-6) O4666 = 6A4P66 — 2p44 - IO5566, '■ 

A longer list is given in Appendix C, while the full set is available online at [15]. 

The same statement is also true for all the 6-index Q-functions, except 
Oii4466- Explicit relations have been calculated down to weight —30. The first 
few are given below, with all available relations online at [15]. 



(-6) 


0666666 


= 40p44 + 


1505566 — 24p66A4, 


(-7) 


0566666 


= 20p36 - 


4p56A4, 


(-8) 


0556666 


= 24p26 - 


12p35 — 204556, 


(-8) 


0466666 


= -20P35 


+ 5 4556 + 16 p46A4 - 2OP55A4 - 8A3 



Proof: Clearly such relations must exist. The explicit PDEs were calculated 
in the construction of the basis, as discussed at the start of this section. 

■ 

Corollary 7.4. There are a set of PDEs that express ^-index p-functions, using 
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Abelian functions of order at most 2. The full set can be found online at [15]. 

(-4) P6666 = 6p^g - 3p55 + 4p46, (28) 

(-5) P5666 = 6p56p66 " 2p45, (29) 

(-6) P4666 = 6p46p66 + 6A4P66 - 2p44 - |p5566 + 3p66p55 + Gpfg, 

(-V P4566 = 2p45p66 + 4p46p56 + 2A4P56 + 2p36, 

(-V P5556 = 6p55p56 + 4A4p56 + 4p36, 



Proof: Apply equation (20) to the first set of relations in Lemma 7.3. 

■ 

The set of equations in Corollary 7.4 is of particular interest beeause it gives 
a generalisation of equation (2), from the elliptie case. A similar generalisation 
for equation (3) would be a set of equations that express the 3- index p- functions, 
using Abelian function of order at most 3. So far the following relations have 
been derived (see [15] for the latest list). 

(-6) pIgg = 4pi6 - '^P56 + 4p46p66 " 8p55p66 " 4p66A4 + 4p44 + 2p5566, 

(-7) P566P666 = 4p|6p56 + 2p46p56 " P55p56 " 2p45p66 + 2p36, 

(-8) p556p666 = -4p26 " 2p35 - 4P55A4 - 4A3 + 2p4,556 - 6p45p56 - 2p46p55 

+ P5566P66 - 2pi6p66 " 2pi5, 
(-8) P566 = 4p56p66 + 4p46p55 + ph + 4p55A4 + 4,p45p56 + 8p26 

+ 4A3 - 2p4556, 

(-8) p466p666 = 4pi6p66 + 4p46pi6 + 2p55pi6 + 4pi6A4 + 2p|g - 4p46p55 

- 2p5,5A4 - 3p45p56 - 2p44p66 - p5566p66 " 2p26 " 2p35 

— 2 A3 + P4556, 

(-9) P556P566 = -2pi6 " 2p45p55 - IP45A4 + p5566p56 " |p34 + |p4555, 
(-9) P555P666 = -3p5566p56 " 4p44p56 + 8p56p66A4 + 12p56p55p66 + lOpfe 

+ f P45A4 - |p34 - |p4555 + 4p45p46 + 8p36p66, 
(-9) p466p566 = -2p25 + f p34 - ^p45>^4 - 2p45p46 + |p4555 + 2p56p55p66 

+ 4p56p46p66 + 4pfg - p5566p56 " 2p44p56 + 4p56p66A4 - P45p55, 
(-9) p456p666 = -P56P55P66 + 2p45p|g + 2p56p46p66 + 2p36p66 " 2p56 

+ 2p44p56 + IP45A4 + 5P5566P56 " |p34 " ^P4555 + 2p45p46- 



Proposition 7.5. There are a set of relations that are bi-linear in the 2-index 
and 3-index p-functions. (See [15] for full list.) There is no analog in the 
elliptic case, although similar relations have been derived in the hyperelliptic 
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and trigonal cases. 

(-6) = -p555 + 2p456 + 2p566p66 " 2p56p666, (30) 
(-7) = -2p446 + 2p455 - 2p466p66 + 2p666'^4 + 2p46p666 

- 2p556p66 + P55P666 + P566p56, 

(^-5^ = -2p56p466 + 2p46p566 + P555p66 " 2p55p566 + P556p56 " 2p366, 
(-8) = -p445 + 2p456p66 + P56p466 - P366 " P566>^4 - P45p666 " 2p46p566, 
(-9) = -2p455p66 + 4p266 + 2p45p566 + 2p466p55 - 2p46p556 + P556p55 

— P555p56 — 2p356. 

Proof: These can be calculated by cross differentiating suitable pairs of equa- 
tions from Corollary 7.4. For example, equation (28) expresses peeee while 
equation (29) expresses pseee- If we substitute for these equations into 

d d 

"5 P6666 — -t: P5666 — t), 

then we find equation (30). 

■ 

A topic of future work in this area would be the construction of relations be- 
tween the p-functions in covariant form, as was recently acheivcd in [9] for the 
hypereUiptic case. 



8 Solution to the KP-equation 

We now demonstrate how such Abelian functions can give a solution to the 
KP-equation. Differentiate equation (28) twice with respect to uq to obtain 

P666666 = 12g|^(p66p666) — 3p5566 + 4p4666- 

Let Uq = x,U5 = y,U4 = t and W{x, y, t) = pee{u). We then rearrange to give 

[W,^^ - 12WW^ - iWt] ^ + 3Wyy = 0, 

a parametrised form of the KP-cquation. In fact, this is just a special case of the 
following general result for Abelian functions associated with algebraic curves. 

Theorem 8.1. Let E he an (n, s)-curve with genus g as given by equation (6). 
Define the multivariate a-function associated with E as normal. Define the 
Abelian functions from a{u) as in Section 4 (with the indicies now running to 
g instead of 6). Finally, define the function W{u) = pgg{u), which we denote 
W{x, y, t) after applying the substitutions Ug = x, Ug-i = y, Ug-2 = t. 

Then, if n> A the function W{x, y, t) will satisfy the following parametrised 
version of the KP-equation. 

{W,,,-12WW,-bWt)^-aWyy = 0, (31) 

for some constants a, b. 
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Proof: Recall Definitions 2.1 and 2.2 which gave the Sato weights for C. These 
can be calculated for the general curve E similarly as 

wt(a;) = —71, wt(y) = — s, 

Wt{Ug)=UJl, Wt{Ug^l) = UJ2, ... Wt(Ul)=LJg, 

wt(Ao) = —ns, wt(Ai) = —n{s — 1), ... wt(As_i) = —n. 

Here n, s are the integers generating the curve, and {wi, . . . ,ujg} is the Weier- 
strass gap sequence for n, s. These are the natural numbers not representable 
in the form an + bs where a, 5 € N. (See [6] Section 1 for more details.) 

Since s > n > 4 we know that {1,2,3} cannot be represented in this form. 
Therefore, wc have 

Wt(Ug) = +1, Wt(Ug_l) = +2, Wt(Ug_2) = +3. 

By equation (26), this implies the p-functions will have weights 

Wt(pff,g) = -2, Wt(pg,<,_l) -3, Wt(pg_l,g_l) -4, Wt(pg_2,g) = -4, 

with all the other 2-index p-functions having a lower weight. Next consider 
Qgggg, whlch wlU have weight -4. This will belong to, T{j,0{2e^a-^'^)), the 
space of Abelian functions defined upon the Jacobian of E, which have poles of 
at most order 2 on O'^"^'. We can therefore express Qgggg as 

Qgggg apg-i.g-i + bpg-2,g, (a, b constants), 

since these are the only Abelian functions of weight —4. We use Remark 4.6 to 
substitute for Q, and then differentiate twice with respect to Ug to give 

Pgggggg = (p^gpg^g) + apg_i,g_i,g,g + 6pg_2,g,g,g. 

CUg 

Then make the substitutions suggested in the theorem to obtain equation (31). 



Further research into the possible applications of these results is currently 
being conducted. 

9 Two-term addition formula 

Theorem 9.1. The functions associated with (7) satisfy the following two-term 
addition formula: 

a{u + v)a{u-v) f, f, ^ 

ayuYayvY 

where f{u,v) is a polynomial of Abelian functions, given in Appendix D 
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Proof: We seek to express the following ratio of sigma functions (labelled 
LHS), using a sum of Abelian functions. 



TTT^/ N <t(u + v)<t(u — v) .„„. 
LHS M, v) = ^ t d ) ^2 ■ 32 

First recall that (7{u) is an odd function, with respect to the change of variables 
It 1-^ [~l]u. We use this to consider the effect of (u, v) i-^ {v, u) on LHS. 

a{v + u)a{v-u) aju + v)a{[-l]{u - v)) 
LHS = / ^2 f \2 = ^ \2 t \2 = -LHS . 

So LHS is antisymmetric, or odd with respect to {u,v) t-^ {v,u). 

Next, recall that cf{u) has zeros of order 1 along 8^^' and no zeros anywhere 
else. This implies that LHS has poles of order 2 along 

(ei^i X J) u (ei^i X J) 

but nowhere else. Together, this implies that we can express LHS as 

LHS = J2Aj{Xji^)Yj{v) - X,{v)Y,{u)), (33) 

j 

where the Aj are constant coefficients (which may be functions of A), and the 
Xj and Yj are functions chosen from the basis in Theorem 7.2. 

Finally we use the fact that sigma has weight +15 to determine that the 
weight of LHS is —30. Hence we need only consider those terms in equation 
(33) that give the correct overall weight. 

We use Maple to construct equation (33) with the Aj undetermined. This 
contained 1348 terms (647 undetermined coefhcients since it is antisymmetric). 
The coefficients were determined using the cr-function expansion. 

■ 

We believe this to be the first of a family of similar addition formula, related 
to the invariance expressed in equation (8). There has been much work con- 
ducted into these addition formula for the trigonal cases (see [10] for example). 
In [14] we see that this has inspired new results in the lower genus cases. 
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A Expansion of the Kleinian formula 

We consider equation (22) from Theorem 5.1. We expand this as one of the 
tends to infinity, to obtain a series expansion in terms of the local parameter 
^. It follows that each coefficient with respect to ^ must be zero, giving us an 
infinite sequence of equations, starting with those below. 

= -Z^ + p46Z^ + (p56W + P26) Z + peew'^ + p36W + pi6 (34) 
= (p45 - P466 - 2w) + ((p55 " P566) W + p25 " P266) Z + (p56 - peee) 

+ (P35 - P366) W + pl5 - Pl66 (35) 
= (p44 - f p456 + |p4666) z"^ + { - Sw"^ + (p45 - |p556 + ip5666) W - |p256 

+ p24 + 5P2666)2 + (5P6666 - ip566 + P46) w'^ + (5P3666 + P34 

— ip356)w + pi4 + 5P1666 - ipl56 (36) 
= (p4566 — |p446 " ^P46666 " 5P455) Z^ + ((p5566 — |p456 — ^P56666 

— ^9555) W - |p246 - 5P255 + P2566 " lp26666)z ~ 4:11)^ + (p5666 " 5p556 

— |p66666 — ^P466)w^ + ( — ^P36666 + P3566 " |p346 " 5P355) W 

— ^Pl6666 + Pl566 " |pl46 " 5P155 (37) 
= -3z^ + (-2p46 - IA4) Z^ + (|p4556 - |p445 " ^P4.5666 " 2p56W - 2p26 

— 3A3 + |p4466 + ^p46^4 + 5^p466666)2^ + ( ^ '2p66W^ + (^p566666 

— f P455 + |p4566 ~ ^P55666 + 5p56-^4 + |p5556 " '^P36)w — ^p25666 

— |p245 ~ 2A2 + |p2466 + |p2556 + ^P266666 — 2pi6 + ^p26^4)z 

(1\5 1 5 5 5\2 

2P66A4 — Y2P56666 + 24P666666 — gp456 + gp4666 + gPSSeej W 

+ (^P366666 + |p3466 + |p3556 + 5p36'^4 — ^P35666 " f P345) W 

+ |pl556 — Ai — ^Pl5666 + ^Pl66666 + |pl466 ~ |pl45 + 5P16A4 (38) 
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B The <T-function expansion 

Wc defined the cr-function expansion as the infinite sum of finite polynomials 

ct(u) ^ Cis + Ci9 + C23 + C27 + C31 + C35 + ■■■ 

We know that C15 was equal to the Schur-Weierstrass polynomial as given in 
equation (25). The other polynomials were calculated in turn using the method 
described in Section 6. The next polynomials, C19 and C23, are given below, 
while the rest of the expansion can be found at [15]. 

Ci9 = A4 • [ 13970880 "6^'^'"4 + xfg Ue^U^^ - Ue'^U^Ui^Ua - U6^M4U2^ + 3^ U5^U2 

+ Jo U^'^Ui^ + ^ Uq^U2^ + 20956320 '"6^^W5^ - \ U^'^Ui - ^ U^^Ui + \ Ui'^U2 

_ 1 „, 7„, 6 ]_ „. 7 4 I 1 13,, 2 _ 1 11 4 , 1 12 

1890 "6 "5 120 6 "4 -r gg5280 "6 "4 ggigg "6 "5 "T 997920 "6 "2 

- gig U6'^W3^ - 5^ Ue^°U4^ + 3^ Ug^U4U2 + \ Ue^Ui^Ui + M6^U5^M4^W2 
" 2^0 "6^W5W4U3 + ^ U6^U5^U4M2 " | U6'*tt5^U4'"3 - f lt6M5^W4^'«3 

- jiy UQ^U^Ui^U^, + i UQUz'^UiU2 - UqU^'U^' + ^ U6'*U4^ + 2^ U6^'"5''w4 
+ ^ M6'^M4^M2 + \ Ue^UiUa'^ + I U5^U3U2 - gfg Uq'^U^^Us + ^ U6^U5''U2 

- M6'^M5^W4^ - 997920 U^^'^U^'^Ui + ^ Uq^Uz'^U2 + ^ ^6^4^ + M6*Ml 

- i U6^U4^U2 - 83Y60 "6"m5W3 " ^ UfC'u^'^Ui^ + ^ U6^U5^W4^ + ^ U5^M4M3] • 



C23(m) = A3 • [^52^32^2 - gig UQ^Ui^Ui + i U5^U4U32 - ^ UQ^U^'^U^i^ 
~ 252 ■"6'^W5^W3^ + ^ U6'*W5"*Ui + i U(iU^'^U2^ + M6^W5^"2^ 
+ Ue'^Uz^U2 + 2^ W6^M5^'"1t|4 M6''w4''m2 " 1796256O "6^^M5''w4 
+ 111767040 Uq^^U^'^U^ + ^;^g^gQ4Q Ue^^U4U2 - ^ U,5^M4U3 " ^ UQ^UzU4Uy,U2 
+ Yg M6'^W4 "1 + 5 "5 U3W2 - 74844 "6 "5 '"2 " g U6 "4 "2 
+ 6M ■"6M5^"4^ + 7^ U6''U5^U4 - ^gig^j UQ^U^^Ui - ^ UQU-^^Ui^ 
+ 2io U6^W5^U4^ + 2T6 W6^W5'*W4'* + Ug'^U^'^Ui^ + U^Ui^Us 

+ 6^ U6^M5''W4^ + 36588 ■"6^°"5^?i4^ + 299376 '^6"w5^W3 " 5 W62u2^ 
- 52o ■"6^"5'^'"4^ + 25147584 Ue^^U5U3 " "6^'"5^W3 + ^ W4^W2^ 

+ 335301120 "6"m4' " g^o "5'"2 - W5'"w4 + ^ U6'"4'u32 

~ T0O8O ^*6 "4W3 + go U6U4°U2 - isfiio "6 ^^4 "2 " igfo ^^6 "4 ^2 
+ ^ Ua^Ui.'^Us + Uz^Ui'^U2 + jfg W6'^"5*^'i2 - 3^ U6^"5''"2 
+ 158760 "6^M5^ " \ UeU5^U4Ui + jg^^ UQ^^Us'^Ui'^ + U^'^Ua,'' + ... 
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' 18144 + 11467298304 ^6 ^"1*5 + 2395008 "6 ''^l + 3 UQU^^UiU^U2 



2 



43 12 11392 1 20 

~ 5987520 "6 '^5'"4U3 + 90720 "6 "5 "4^2 — 229345966080 "6 ^4 

- ^ U6^M5^U4^'«3 + 5J0 ■W6^U5*'"4^ + 5 U6^M5^M4M2^ + lygio W6''"5''u4 

1_ „ 8,, 5 7_ ,, 4,, 5 _ 1 ,, 11„, 4 _ 1 „, 15,, 4 

8064 "6 "4 360 "5 "4 3991680 6 "4 75442752 "6 "5 

1 „, 11„, 2 I 1 1, 16„ 4_ J_ „ 5,, 6 1 „, 14,, 3 

199584 "6 "3 28740096 6 "2 720 "6 "4 167650560 6 "4 

- U6^M5^U4U32 - Uq^UzU^U2 + 5 U6^U5^"4^Wl + M6^W4Wl 

+ I U6^W5W4^U3 + U6^U5^U4W3 + ^ U6^W52u4Ui - ^ U(iU^^UiU2 

- gig U6^U5^1i4U2 - 3^0 Uq^Uz'^U4^U2 ~ Uq'^ Ur^U^^ U3 - ^ Ug'^Ui 

- ^ U6''m5^U4W3 - T2 "6^U5'*W4^W2 + 45^60 '«6^'"5tt4^W3 + 5 U4'*Ml 
10595783632896 "6 ^ " g ^6 W4'^W3 + g U6^ti5 U4^U2 ~ 641520 "6 ^5^] 



' -^4 • [ ^ TM "6''U4^W1 - 7^ Ue^U4'^U2 - 204^200 U6^^W5'^?^4 
- 2993760 '"6"U4'"2 + 3Y5 U^'' + 22^ U6^U4Ml + 5 UQ^h 
' ^ U5^U3U2 + 467^775 W6^^W5^W2 " 5 W6^W4^'«2^ - ^ M6W5' 



. + i U6^U4^Ui 

i-^D '-^^^"z I 467775 ""D "'O 3 "'O "'4 "'^ 315 ^^^^ ^4 

■ y|g U6''U5^U4 - ]3J75 U6^U5'5u4 + u^U^'^Ui^ + Ufi^U^'^Ui^ 

' U6^U5^U4'' - 2I0 UQ^U^'^Ui'^ - ^ U^Ui^Uz - W6''W5^W4^ 

I 101 „, 10„, 2,, 3 4 11„, 3„, 8 „, 2,, 6„, 3 

+ 226800 "6 M5 "4 - 467775 "6 "5 "3-135 ^6 M5 -U4 

- 4^ W6^M5^?^3 + I U4^W2' + 83825^2800 "6'^M4' + M5*W2 

+ ^ Ue^Ui^Us'^ ~ 2^20 "6*^W4U3^ + M6U4^"2 + jg^oo W6^°W4^"2 
+ sfio UQ'ui^U2 - gig U6^W5^M3 + i 1152^4^112 + jig U6'*?l5^"2 

1 2„ 3,, 2 43 7,, 8 1 13,, 2,, 2 , 1 2,, 7 

- 3 U6 U4 U3 - 396900 "6 "5 " 199584 "6 "5 1^4 +34 "6 "4 

26 192 1 19 2 1 12 

- 45 W5 Ul — 22680 "6 "2 — 28668245760 "6 "5 " 14968800 ^6 "l 

+ W75 "e^Wsl" - I UQ^U^'^Ui^U^ - ^ U6U5^W42ji3 - 74^^ W6^^"5W4W3 

+ 2^5 "6^M5'?i4W2 + 573364915200 "6'°"4 " ^ U6^M5^?^4'm3 

1 „, 5„ 6,, 2 , 29 „, 9,, 4,, 2 247 „, 8,, 5 , 1 4,, 5 
~ T350 "6 "5 + 68040 "6 "5 "4 — T0O8OO "6 "4 +35 "5 "4 

593 ,, 11„ 4 I 1 1, 15„ 4 I 1 11„ 2 I 1 ,, 16„ 

9979200 6 "4 T 67359600 6 "5 1247400 6 "3 2514758400 6 "2 

1."^ ^fi. 1Q 1/1'? 1 /I 



593 ,, 11„ 4 I 1 1, 15„ 4 I 1 11„ 2 

9979200 6 "4 T 67359600 6 "5 1247400 6 "3 

I 13 „, 5„ 6 I 19 „, 14„ 3 1 ,, 4„ 
+ 900 "6 "4 + 41912640 "6 W4 - 13 "4 "i 

+ M 'U-&^UiU2^ - ^ ?l6^'«5U4'*W3 - 9I5 U6^U5^U4U3 + ^ ?l6U5'^M4U2 

+ ^ M6^M5^M4U2 + ^ Uq'^ U^'^ Ui'^ U2 - 3^ U6^"5W4^M3 + | UQ^U^'^Ui^h 

- ^ M6'''"5^M4U3 + 3 W62m5''?^4^'«2 - 2^0 UQ^U^Ui'^Us 

1 „ 11,, 6 1 ,, 23 I 1_ „ 8„ 4 

561330 "6 "5 26489459082240 "6 "T 1890 "6 "5 "2 

39916800 ■"6^^"5^U4 — 15717240O "6^^W5?^3 " 4725 ^5"^°"4] • 
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C The 4-index Q-functions 



This appendix contains a list of PDEs, expressing the 4-index Q-functions that 
were not elements of the basis for r(j, O(20[^l)), as a linear combination of 
basis elements. This appendix contains all the equations down to weight —22. 
The full set can be accessed at [15]. 

The PDEs are ordered in decreasing weight (as indicated by the number in 
brackets) . 



(-4) 
(-5) 
(-6) 
(-7) 
(-7) 
(-8) 
(-8) 
(-9) 
(-9) 
10) 
10) 

■11) 

11) 
.12) 
.12) 
.12) 

13) 
13) 

13) 
.14) 
.14) 



Q6666 
<35666 
Q4666 
Q4566 
Q5556 
Q4466 
(35555 
Q4456 
<33666 
Q2666 
Q4446 

Q3466 
Q4445 
Q2466 
Q3555 



V2555 
(33446 

(33455 
(3l666 
(32446 



(-14) Q2455 



-3p55 + 
-2p45, 

6A4P66 - 
2A4P56 -f 
4A4p56 ^ 
4A4p46 ^ 
I6A4P55 
|A4p45 - 
2A4P45 



4p^ 



5566) 



46, 

2p44 - f ( 
2p36, 
4p36, 

A4P55 + P35 
f 4p35 + 24p26 - 6Q 
- 2p25 - |p34 

^(345557 



6p26 - (34556 + 4A3 



4556 



4555 J 



?4455 — 5(33566 — 5A4(35566 + ' 
,2, 



I6A3 



366A3 



6A4P44 - 2p24 - 2A4P66 

— (34455 + 1-^4(35566 ^ f (33566 + 4p66A3, 
4p36'^4 — 5(335667 

6p36A4 — 2p56A4 + 8p56A3 — 3(32566 — §(33556: 

8pl6 - P33 + 2A2 - 5(32456 - (33456 + 4A4p26 + 2A4P35, 

24pi6 - 8p33 + I6A4P35 - 6(33456, 

-12pi6 + 9p33 + 6A2 + 12Q3456 + 12p55A3 + 4p46A4 

— 16p46A3 - 3P55A4 4- I2A4P26 - I8A4P35 - 6(32556, 
2pi5 - 8p23 + I6P25A4 - 6(32456, 

6pl5 - 4p23 + 4P25A4 + IP34A4 
+ |p45'^4 ~ 5(34555-^4 — 2(32456, 

— 4P25A4 -I- IP34A4 — IP45A4 + ^(34555.^4 + 2(32456, 

— h^4Q3566 + (33366 ^ 5 (33445, 



6 

|Pl4 



2p22 + IP24A4 



3p66^4A3 



4455 — 
lp66'^4A3 



^(35566A4 -|- 2P44A3 

2p22 + f P24A4 — |p66'^4A3 — §(33366 
^(35566-^4 + 6P66A2 — 2P44A3 + ^A4 (34455 
^(33445, 



qQsSGG + (35566-^3 
A4 (33566 — |(33445, 
(35566-^3 
-^4(33566 
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(-15) Qi,566 — —jQ2445 ~ 1(52455^4 + |p2366 " |p23p66 

3 , 



3p26p36 



■15) Q3356 

•15) (53444 

■16) Q3346 

■16) (33355 

■16) (32444 



■17) (3l456 

■17) Q1555 

■17) O2355 

■17) O3345 

■18) Ql446 

■18) (32246 

■18) O2255 



■19) (3l366 



■19) Q3336 



— 5p36-^3 + |p56A2, 



40 



2445 



SpaeAa + 3p56A2 — -^(32566-^4 — ^(3 



— jQ2445 + 4p36A4 — §(33556^4 + ^p36A3 



2^2366: 
|p56A2 



^2566^4 — x02366, 



— 2(3l466 + 4Ai — (3l556 ^ (32266 + I6P16A4 + SpasAa — P55A2 
+ 2p26A3 — 2(32356, 

8(3l466 - 8A1 + 6(3l556 - 32P16A4 + 4P35A3 

— 2(32356, 

+ 2A4A2 + lOAi + 36P16A4 + f P35A3 + ^P55A2 - (33456A4 

- |(32556A4 - 7p26A3 - 6p46A2 + 3P33A4 + 6A4P3456 

- 12A4p34p56 - 12A4p35p46 - 12A4p36p45 + 4:p26^l " 



'35 A4 



1466 — 9(3l556 + 2'32266 + 3(32356, 

- -2pi3 - P45A2 + 2p25p26 + P22p56 + 4P15A4 - (32346 

— tl'54555A3 — 5P2256 + P25A3 + 5P45A4A3 + IP34A3, 

- 4pl3 + 6P45A2 + 3(32256 - 8P15A4 + 6(32346 + |'34555A3 

— 6P25A3 - 2P45A4A3 - 8P34A3, 

: -4pi3 - 4P45A2 - 2(32256 + 8P15A4 - 2(32346 
+ 4P25A3 + 4P34A3, 
: — 2P25A3 + IP34A3 — IP45A4A3 + 2P45A2 + (32256 

- -6pi2 + 6p66Ai + |(33566A| - Qu55 + 8P14A4 - 

— (32345 ^ 5Q3344 ~ 5'33366A4 + ^Qsiid^i^ 
■ 2P44A2 — 4p66Ai — 2P66A4A2 + 5(33566A3 + Q3366^4 

tA4 — (33566A4 + 2P24A3 + |(35566A2 + (3l455, 



- ^'34555A3, 
3566 A3 



^3445^^4 



- -8pl2 - 8P44A2 + 16p66Ai + 4^pi4A4 + 8P66A4A2 - 2p2345 
^ IQ3566A3 — |(33366A4 + |(33445A4 + ^Q3566^4 + |(34455A3 
+ 3A4A3P5566 — IA4A3P55P66 — |A4A3p|g — IP66A3 + IP24A3 



+ 4p23p45 + 4p24p35 + 4p25p34 - 
t(32245 ^ 1(32344 " -j7T(32566A3 



3(35566 A2 — 2(3l455, 



+ m'32566A| + ^p56Al — ^p56A4A2 



P36A2 — i:(33556A3 



_3 
' 10 
_3 
20 

9 

5 

-to 

_3 

10 

#P36A4A3 



f2366^ 



A4 + ^(32455A4 



T7TP36A4A3 



2245 + f (32344 — |(32566A3 + |( 
17l(32566A4 + ^p36A2 — ^p56Al - 



b366A4 + ^(32455A4 

3 
5 



P56A4A2 



;3556^3, 
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(-19) Qi445 

(-20) Q2244 



hQ 



2245 

^(32566^4 I 2^ 



(52344 — 5'52566'^3 " §(32366^4 + SpseAi + 3p36A2 
5*52455^4 — 2P56A4A2 + 3P36A4A3 — IQ3556A3, 
+ 2A3A2 + 28A4A1 + 2P55A4A2 - 6P35A4A3 + 3P33A3 

— 24p46Ai + 112pi6A4^ - 44pi6A3 + 28P55A1 + 9P35A2 

+ 2OQ1356 — 6(5l444 + IOQ1266 — 22(5l466A4 + Q2266A4 
+ 2(52356 A4 + 2(53456 A3 — 22(5i556A4 — (52556 A3 — 6P26A2, 
8PI6A3 — 2p55Ai + 2P35A2 — 2(5i3,56 — 2(5i266, 

-32A4A1 + 24p46Ai - 128pi6A4^ + 56P16A3 - 32p55Ai 

— 30(5i356 + 8(5l444 ~ 12(5l266 + "^^Qli&e^i + "^^Qlbb^^i 
+ 4P35A2, 

(-21) (5l256 = ^5(52236 ^ (5l346 ^ ]^(54555A2 + ^p45A4A2 + 3P15A3 
+ IP34A2 - 2p45Ai, 
\Q22m + (5l346 + g:(54555A2 — \ 

— §(52236 — 3(5i346 + \Q4:555^2 
+ I2P34A2 - 8P45A1, 



(-20) Q2336 
(-20) Q3335 



(-21) Q1355 
(-21) Q3334 



i;2335 
§(52335 



P45A4A2 + P15A3, 

P45A4A2 + 9P15A3 



(-22) (5i345 = SpeeAo - pn + jf pi4A| + ^p66A4Ai gg 



I (53344 A4 
?2345A4 



^(53566 A4 

2 

5 



2-^1255 

3 J_ 

15 



^Q3566A4A3 
^(5l255 + ^Pl4A3 — P44A1 — (5l246 

3366A4^ — \Qd3QQ^d + ^Qmb^i 
P12A4 + ^(53445 A4^ 



(-22) (52226 — — 6pii + §(53566A2 + \QAi55^2 + 6(55566Al 



-P14A4 



-^p66A4Ai — f§(53566A4A3 — ^pl2A4 + ip24^2 — lOpeaAo 



I2P14A3 4 

3344 A4 
(53566^ 



6P44A1 



6 
5 

4/n \ 3 



2P66A3A2 

2 , 3 



(53366 A4 + 2^ 
^ (55566 A4A2, 



6(5l246 
3366A3 — |(^ 



^ (53445 A4^ 



' 5 
'3445 



A3 — '^(52345A4 



(-22) Q 



2235 



-2pii — (53566A2 — 4(55566Al — fipl4A4 + ^p66A4Ai 



f|(53566A4A3 + IP12A4 + 10p66Ao — Q 



1255 -t- IP14A3 
-^(53445A4^ + i(53344 A4 + jrQsSGe^ 



IOP44A1 
' ' =;A3 



2Q 



1246 



^QssGG^a + ■gQsiib^s + 5(52345 A4 " 



15' 

15^3566^^4 , 
32 , 



(-22) (52334 — — 4pii — (53566A2 + 2(55566Al — Y5P14A4 5-P66A4A1 



£ 

15 



'3566A4A3 + IP12A4 + (5l255 + IP14A3 + 4p44Ai 



(53445A4 + i(53344A 



^( 

15 



^3445 A3 



A3 + 



|(52345A4 — 1%(53566A4 



3366 A4^ 
3 



§(53366 A3 
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D The two-term addition formula 



The Abelian functions associated with C satisfy the following two-term addition 
formula. 

a{u + v)a{u-v) n n 

where f{u, v) is a finite polynomial of Abelian functions. We write f{u, v) as 
/(m, v) = P30 + P26 + P22 + P18 + Pii + Pio + P6 + P2 

where each Pk contains the terms with weight — fc in the Abelian functions and 
weight A: — 30 in A. 

P30 = jQll4466{u) + lQ3566{v)Ql356{u) - | pl4 (■u)(32356 (") 
+ ^Ql556{v)pl4{u) - i(52236(M)p25(^^) " |Q2345(m)P33(^^) 
+ ]^Q3344('y)Q3456(M) " 5 P25 (■y)(3l346 (") + 5^4556 (^')pll (") 
+ |(93344(M)pl6('y) + JlQl346{v)Q4555{u) + ^Qll45{v)p56{u) 
+ j2Q2236{v)p34{u) + ^Ql444{v)Q3566{u) ~ P46 (^^)Qi146(m) 

- |Qi556(«)P22(m) + |(3ll55(w)p46(^') + |Q3566 (l')(5l266 (") 

+ |P12('w)Q3456(m) + P26(«)Qi246(m) + Q2346 (■w)pl5 (m) - 5P12(^')(32556(m) 

+ ^02345 (M)Q2556('y) " |Qi356(?^)(34455 (") ~ ^ Q 2245 {v)Q 2566 {u) 

+ P36{v)Ql245{u) - |p33(^')Q3344(m) + Ql246 (■w)p35 (") " |P24(^')Qi444(m) 

+ |<3i466(m)Q3445('w) " |(3i466(m)P14('w) " P26(«)Pii(m) " (32346 (■y)p23(M) 

+ |pl6('y)Q2345(M) - ]^Q3445('y)Q2266(M) + ^Q5566 (^^)Ql244 («) 

+ ^Ql466{v)Q3366{u) + j8Q4555{v)Q2335{u) + 3Q44.55 (■w)Ql266 (") 

- ^(34455 («)(3i444(m) + |Q2345 (■w)'93456 (") + 5Q1156 (■y)p45 (m) 

+ Q2256('y)pl5(M) - ^P12(^^)P16(m) ^ lQl356{v)p24{u) + ^Q2335 (^^)p25 («) 

- 5p33('y)'9l455(M) + 5P44('w)0i166(m) " P35 (■w)pll (m) - Qi455(«)P16(m) 

- ^(52236 (■w)(34555(m) " ^Q3344('«)Q2556(m) " ]1j(32566('w)Q2344(m) 

- ]^Q3556('y)'92245(M) " ^Q3556('w)'32344(m) + |Q2445('^)Q2366(t*) 
+ ^Ql556{v)Q3366{u) - \Q\255{u)pi^(v) - \QlX44{v)p66{u) 

- ^Q2335{v)p34{u) + |Q3366('w)(52266(m) + |P24(w)(3i266(m) 

- i(32256(«)P23(m) - P22(m)Qi466('w) + 5 P55 146 (") + P44 (^')(3l244 (") 
+ ^Ql346(M)p34('y) - |Pi4(m)Q2266(^^) + Pl3 (■y)Q2456 (") + ^P33(m)P12(^^) 



/2356 



3v,;z.30d(^')(33366(m) + 5Q3445 (■w)'9l556 (") " 1^3445 (■y)(32356 («)• 
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^26 = [jQll55{u) - (5i146(m) - ^Q25.56 (m)'53366 (■«) " ^Ql556{u)Q4i55{v) 
+ H<33566(w)Q2266(^) - |Q5566 (■w)'9l266 (") + 5(3l466(«)Q4455 (") 

- J^Q3445{v)Q3456{u) - |(95566 (■y)Ql356 (") + IQ 3566 {u)Ql 556 (v) 
+ ^Q2556(w)Q3445(^') - 5Q5566(m)Qi444(i') + ^OsSSe (■w)'52445 («) 

- m'33456('w)Q3366(m) + ■^Q2566{u)Q2A45{v) + |Qi466 (■w)'53566 («) 
+ ^Q3566(w)Q2356('y) " Pll (m)P55 (^') " ^Q3445('y)p33(M) 

+ ^Q2556{u)pl4:{v) - ^p45(M)<3l346('y) " ^P16('w)P14(m) 
+ 2p23{u)pi5{v) - ^Q3366(«)P16(m) ^ ^Q2335 (m)P45 (■«) 

+ |P33(m)(33366(^^) + 2(5i466(^^)P24(m) + 4pi6 (f )p22 (") ^ 2gi556(M)p24(f ) 
+ 3P44(m)Qi266(^^) + |p44('y)(9l444(M) + |p36(«)Q2344(m) 

- 2p35(M)pl2('y) + ^p45{u)Q2236{v) - J^Q3456{u)pl4{v) 

- p66{v)Ql244{,u) - 2p25{u)pi3{v) ~ |p36(m)Q2245(^^) 

+ ifpie WQ3445(m) + 6P44(m)(5i356('w) " 33(M)pl4('y)] A4 

-P22 = [p35('y)pl4(M) - Ql466(v)Q5566{u) - P23 )p25 (") ^ f|pl6 (m)(93566 (^') 

- Ql556{u)Q5566{v) - ^(32566 (m)Q3556 (■") " ^ Q 3456 {u)Q 3566 (v) 

- 2Qi466(m)P44(^^) - 5Q3366('y)p35(M) + 2p44(M)Ql556 (^^) 
+ h^l5{v)p25{u) - J^p23{u)Q4555{v) + I p33 {u)Q 3566 {v) 
+ ^Ql266{u)p66{v) + ^Ql356{v)p66{u) ~ |p23(^^)P34(m) 

+ ^Q2556(m)'93566('w) + \Q\444{u)p&(,(y) + ^(33445 (^')p35 («) 

- |pl6(M)Q4455('y) - 4pi6(M)p24(^^) - iQ2366(w)p36(^^) " T ^^15 (■")p34 (f) 
+ 5iPl5(M)Q4555(^^)] A3 + [ - ]^Q3456(m)Q3566(«) " 4p44(M)(3l556('y) 

- |(3i466('w)(55566(m) + 3Qi266(m)P66('^) + \Q2366iu)p^Q{v) 

- l^Q2445(w)p36(«) +4Qi466(M)p44('y) " ^P33(m)Q3566(«) 
+ ^Q2556(w)Q3566('y) " |Ql444(M)p66('y) + ^Pl6(M)p24('y) 

+ |P16(m)(94455('w) + ^pl6(M)(33566('y) + ^(32566 (m)(33556 (l') 

- 6Qi356('y)p66('") + iQl556(M)(35566(^')] A4 
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Pl8 = [juQ2566{u)p36{v) - p23{u)p45{v) - ^Q2345{v) - JqQ3344{v) 
+ |P12(m) + |p56(M)(52366('y) + ^'34555Np25(^') + |P34(m)P25 (^') 

- ^p34(v)(34555(«) + ^P55{u)pl4{v) ~ 2p4r,{u)pir,{v) - ^p55 (m)Q3445 ('^) 
+ |P55(m)Q3366('^) - fQ3556(M)p36(^')]A2 + (52566 (w)p36(w) 

- f Q3556(m)P36(i') + 2p66(M)Ql556(f ) + 2p66(M)Ql466(«) + |p44('w)pl6N 

- ipl6('w)Q5566(u) - ^p45{u)pl5{v) + ^ p23 (w)p45 (f) 

- iQ3566(M)p35(^')]A4A3 + [4p66(m)Qi556('w) + ]^Q2566(m)P36('w) 
+ 8Pi6('u)(95566(m) + 4p66(M)(3l466('y) " f^p44(^^)pl6(M)] A| 

Pl4 [P25(m)P45(^^) + P22(m) + |Q3445(m) + ^Puiv) + J^QsaGsiv) 

- ^Q3556{u)p56{v) ~ 6p44{u)p26{v) + 8^5566 (■u)P26(m) + ^Pi6(m)P66('w) 

- ^<34455(^')P46(m) + f|Q3566('w)p55(M) ^ ^Q4555 (m)P4.5 (^') 

- fip33Np66('w) - |P24(m)P55('w) " f (93456 (m)P66('w) " 2p24(l')p46N 
+ |(94455(«)P55(m) - 3p35(M)(95566(t') - |p34(m)P45(i') " 3-P35 (M)p44('y) 
+ ToQ2556{u)p66{v) - ^Q3566{v)p46{u) + ]^(32566(«)p56(lt)] Ai 

+ [^(93566(^')P55(m) + ^p25(M)p45('y) " 1^(53556 (w)p56 (^') 

+ ^p34(M)p45('y) - ^(92566(^')P56(m) " ^(93366('y) " ]^(93445(m) 

+ ^Pl4(^')] A4A2 - ^Pi6(^')P66(m)A| + ^pl6('y)p66(M)A4A3 

- f Pl6(^')p66(M)Ai 

PlO = [ - 2p24(t') - ^Q3566{v) - | P44(t))p55 (u) + 5p3(i{u)p5(i{v) 

- 54(94455('y) - ip35(«)p66('w) + |(35566(w)p55(l')] Aq + [4p44(t>)p46(M) 
+ |p24(w) + ^P36(m)P56(«) + 6p26('w)p66(M) " ^^p3^{u)pm(v) 

+ 6p44('u)p55(M) - |(95566(W)P46('W) + |(95566(m)P55 (f) " §(93566(m) 

— |(34455('W)] A4A1 

A = ^p66(M)p55('y)A3Ai + ^(55566(^')A4Ao - 4p46(M)p66(u) A|Ai 

- 4p46(v)p66(M)A3Ai + |p55(M)p66(f )A4Ao - 6p55(t>)p66(M)A|Ai 

- §(95566 (■") A3 Ai + |p44(t))A4Ao - ^p44(v)A|Ai - 2Q5566 (f ) A| Ai 
+ 3P44(m)A3Ai 

Pi = |p66(m)A4A3Ai - f P66(m)A2Ai - ^p66(M)AtAi + |p66(m)A1Ao 

- f P66(m)A3Ao 



31 



